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INTRODUCTION 


Models can be formulated at different scales, in particular: 

- the micro-scale of individual particles of the solid phase and of 
individual pores; 

- the meso-scale of individual roots, aggregates, and “hot spots" of 
intense biochemical activity due to the presence of organic matter; 

- the macro-scgle of the entire, potentially rooted zone at a scale of the 
order of 1 m 

- the mega-scale of the field; 

- the giga-scale of the region. 

Soil physicists feel most comfortable facing problems at the 
macro-scale, This scale offers the best perspective to formulate models in 
terms of measurable quantities. Most methods for determining physical and 
physical-chemical properties pertain to the macroscopic scale. Most 
instruments used to monitor actual conditions in the field deliver 
information at the macroscopic scale. 

The macro-scale is the natural basis for integration to the mega-scale, 
i.e, the field scale of concern to the farmer, At the mega-scale attention 
focusses on integral inputs, storages, productions, depletions, and 
outputs. 

The macro-scale itself is based on integration from the micro-scale 
either directly or indirectly via a meso-scale. The growth and acitivity 
of organisms is strongly determined by the physical and physical-chemical 
conditions at the miero- and meso-scales, Therefore there is a real need 
to study conditions and processes at these smaller scales. Two approaches 
should be further developed. One is the development of models and 
experimental techniques appropriate to these smaller scales. The other is 
to explore the potential of inferring from measured conditions at the 
macro-scale the conditions at the smaller scales. 

This paper reviews some aspects of the origin and roles of length and 
time scales associated with movement of water in unsaturated soils. 


THE PHYSICAL MATHEMATICAL MODEL 
The balance of mass for the water may be written as 
a0/8t = - V. (ðv) - u, [1] 
where t is the time, y is the vector differential operator, 0 is the 
volumetric water content, y is the velocity of the water, and u is the 
volumetric rate of uptake, The volumetric flux êy is given by Darcy's law: 
éy = - kih] vh + k{h] vz, [2a] 


v = - D{#] V9 + kle] Vz, [2b] 
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ive - vg + k[¢) Vz, [2c] 
where h is the tensiometer pressure head, z is a vertical coordinate with 
its origin at the soil surface and taken positive downward, and the 
diffusivity D and the matric flux potential ¢ are defined by: 


D = k/C, Ëi fen nf kan a p, 200, [3] 


where C = d@/dh is the differential water capacity, Symbols in square 

brackets denote functional dependence. Unlike the dependence of k upon @, 

the dependence of h upon @ is subject to hysteresis. As a consequence, 

Eqs. [2b] and [2c] are, strictly, only valid for monotonic changes in 

water content from some initial condition with uniform @ and h. 

Two groups of parametric expressions describing the relationships 
between water content, pressure head, and hydraulic conductivity have 
received a lot of attention: 

I A group yielding flow equations that can be solved analytically, in 
most cases as a result of linearization following one or more 
transformations; 

II A group that is favored in numerical studies and to a large extent 
shares flexibility with a rather sound basis in Poiseuillian flow in 
networks of capillaries. 

To group I belong the following classes of soils (Broadbridge and 
— 1988; Raats, 1983, 1988; White and Sully, 1987): 

Linear soils, with h proportional to ln ((9 - 2 WALAN - >) and k 

proportional to (@ - 6@,), where the subscripts I and u denote lower and 

upper limits, respectively. 

2. Green-Ampt or delta — soils, with the diffusivity all 
concentrated at f = @ ; 

3. "Mildly" nonlinear soils, with h a linear function of 9 and k an 
exponential function of @, implying exponential k{h], D[@], and D[h] 
relationships and leading to linearization of Darcy's law in terms of 
the matric flux potential; 

4. "Versatile" nonlinear soils, with D[#@] and k[@] chosen such that the 
one-dimensional Richards’ equation can be linearized by successive 
application of the transformations of Kirchhoff, Storm, and Hopf/Cole. 

Classes 1.3 and 1.4 include as extreme subclasses the classes I.1 and I.2. 

Class 1.3 has been useful in obtaining analytical solutions for 

multidimensional steady flows and for steady flows involving uptake by 

plant roots, The solutions of the Richards’ equation for class 1.4 

represent a major advance in theoretical soil physics since the last ISSS 

congress (Broadbridge and White, 1987, 1988; Broadbridge and Rogers, 1990; 

Sander et al., 1988; Warrick et al, 1990; White and Broadbridge, 1988). 
Group II can be best represented by a superclass of soils based on (see 

Raats 1990a for further details): 

a) A four-parameter relationship between the pressure head h and the 

degree of saturation S = 8/8: 


} <9 Y 
h/h og * 2 J'y i4] 


where h £ is a characteristic pressure head and a, J, and y are empirical 
constan 


b) A three parameter dependence of the relative hydraulic conductivity 
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k/k; upon h and S: 


Sras 
a ET as š 
k/k TS í Tara ) [5] 
0 


where 1, p, and q are constants, Introducing [4] into [5] gives 
k/k, = S7 (Iga ((pAY + 1)/a, (1 - py)? [6] 


where Ia | . , . } is the incomplete beta function. 

The Superclass of soils, in sofar it is defined by [4], is related to 
well-known subclasses (see Raats, 1990a for details). Equation [5] reduces 
to Burdine’s expression if r=n+m - l, p= 2 +b, and q = l and to 
Mualem’s expression if r=n+m - 2, p= 1l +b, and q = 2, where n is the 
connectivity parameter with 1 < n < 2, and m and b are tortuosity factors, 


INTRINSIC CHARACTERISTIC LENGTHS, VELOCITIES AND TIMES 
At any state of an unsaturated soil characterized by particular values 


of 9, h, and k one can define the differential intrinsic characteristic 
length, velocity, and time by 


A = kdh/dk = dé/dk, [7] 
u = dk/dé, [8] 
r = A/v = k(dh/dk)(d@/dk) = (dø/dk)(d9/dk). (9) 


In general these differential parameters will be a function of the water 
content. However, corresponding robust integral parameters can be defined 
by averaging A over k, and v over 6: 


h 

— gy 
ys tk — [10] 
. % =k 
=< (11) 

u ees A h 

(69 - 6), f “keh 
— 1 h (6 - 0.) - $) 
r=- 1 ete [12] 
(kh, > kD) (Ko ED 


The differential parameters À, v, and r were considered earlier in the 
context of "mildly" nonlinear soils, for which A = constant, v ~ k, and r 
~ k ` (Raats, 1976). * 

For h = - œ, k ~ 0 andh = 0, k =k, the negative of à corresponds 
to the chitical préssure head"h_ of Bouwer (1964). Raats and Gardner 
(1971) derived expressions for fi, for six empirical relationships between 
k and h. 

The integral parameters à, v, and r play a key role in numerous recent 
papers from the CSIRO Centre of Environmental Mechanics at Canberra (e,g. 
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White and Sully, 1987, Broadbridge and White, 1988), In those papers A is 
written in the form 


2 
S 


Nap a 
CO, - Ak -ka 


[13] 


where S[h,, h ) is the sorptivity and b is gbtained by comparing 

[10] and 13) "and using an expression for S” based on the quasi-analytical 
solution of the horizontal adsorption problem of Philip and Knight (1974). 
For the resulting expressions for A, 0, and F, associated with the classes 
of soils belonging to Group I, the reader should consult the papers of 
White and Sully, 1987 and Broadbridge and White, 1988. 


THE MILLER SCALING THEORY 


The scaling theory of Miller and Miller (1956) is concerned with 
geometrically similar media characterized by length scales à „= l and A 
For geometrically similar media with geometrically similar distributions 
of water and of air, the STVF (= surface tension, viscous flow) -theory 
implies that simple relationships exist between water contents, pressure 
heads, and hydraulic conductivities: 

1. Geometric similarity implies @ = 6,; 

2. The inverse relationship between the pressure head and the mean radius 
of curvature implies h = A hy; 

3. The linearized Navier-Stokes equation at the microscopic scale implies 
that in Darcy's 8 at the macroscopic the hydraulic conductivity 
satisfies k = \ ee 

The three primary scaling rules just given are of the form (Raats, 
1983): 


wo Am Wy [14] 


with integer n. Using the three primary scaling rules, secondary scaling 
rules can be inferred from Darcy's law and from the balance of mass for 
the water. Darcy's law implies simple scaling rules for the spatial 
coordinates x, y, and z, the velocity y, the volumetric flux fy, the total 
head H = h + z, the diffusivity D, and the matric flux potential ¢. The 
volume balance for the water implies scaling rules for the time t, and the 
volumetric rate of uptake u. Scaling rules for the differential water 
capacity C, characteristic length A, and characteristic speed v, all three 
potentially a function of the water content, can also be inferred easily. 
Table 1 of Raats (1983) gives the values of n in Eq. [14] associated with 
various parameters. Most noteworthy are the scaling rules for the spatial 
coordinates and the time 
- the spatial coordinates, and hence all macroscopic length scales in 
processes, should be inversely proportional to the microscopic length 
scale; 
- the time coordinate, and hence all time scales in processes, should be 
inversely proportional to the cube of the microscopic length scale. 

The various parametric expressions for the relationships between 6, h, 
and k (Groups I and II above) can be interpreted so as to be consistent 
with the scaling theory of Miller and Miller. For example, if equations 
[4] and [5] are so interpreted, then hog is inversely proportional to An 
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and k_ proportional to A_. = 

An important implicatfon of the power function dependence À of all 
the variables w in Table 1 of Raats (1983) is that if the length scale À 
is lognormally distributed, then all the variables w/w, will also be 
lognormally distributed. This is a consequence of the reproductive rule 
for lognormal distributions: if the variable x is lognormally distributed 
with mean u and variance o then e x is lognormally distributed with 
mean a + by and variance (be)”. 


DISCUSSION AND CONCLUSIONS 


The characteristic lengths and times discussed so far are intrinsic 
soil properties. In actual flow problems two other categories of lengths 
and times arize 
- lengths and times associated with geometry and boundary conditions; 

- lengths and times emerging as consequences of a flow process. 

Examples of length scales associated with the geometry of the region in 
which the flow takes place are: the length of a soil column, the depth of 
a water table, the location of a drip irrigation source within a flow 
region and/or relative to other sources. Length scales may also arize from 
boundary conditions: the size of a disc source imposed at the soil 
surface, the characteristic length associated with a pattern of wetting, 
Boundary conditions may also imply characteristic times; the duration of 
an input at a certain rate, the period of a cyclic input. 

Examples of lengths emerging as consequences of a flow process are the 
maximum height of capillary rize of a given upward flux from a water 
table, the size of a pond at the soil surface resulting from an imposed 
flux, the effective thickness of the capillary fringe with regard to 
lateral flow, and the maximum thickness of the zone depleted ahead of a 
moving root front. 

The characteristic times emerging from flow processes are of two types. 
One type is the time elapsed when a process settles down to some 
asymptotic behaviour such as a steady state or a time invariant travelling 
wave, or a regime in which the growth or decay of the profile is a simple, 
explicit function of time. Another type marks moments at which imposed 
flux boundary conditions become untenable. Three important examples of 
this in soil physics are: 

1. The rate of evaporation imposed by the evaporative demand of the 
atmosphere cannot be sustained and switches to a falling rate. The falling 
rate will be inversely proportional to the time lapsed since some properly 
chosen instant. 

2. The rate of infiltration imposed by patterns of rainfall or irrigation 
cannot be sustained and switches to a falling rate. The time at which the 
switch occurs is called the ponding time. This ponding time turns out to 
be a robust parameter, e.g. being nearly independent of soil type within 
the wide class of "Versatile" nonlinear soils. 

3. The rate of uptake of water imposed by plant demand cannot be sustained 
and switches to a falling rate, thus leading to deficiencies. 

In many cases there emerge in the formulations and solutions of flow 
problems ratios of the intrinsic characteristic lengths \ and i and times 
r and 7 on the one hand and the characteristic lengths and times 
associated with geometry and boundary conditions or emerging as 
consequences of a flow process, 
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SUMMARY 


This paper reviews some aspects of the origin and role of length and 
time scales associated with the movement of water in unsaturated soils. A 
distinction is made between differential and integral scales. The links 
with parametric expressions describing the relationship between water 
content, pressure head, and hydraulic conductivity and with the Miller 
scaling theory are indicated. Also discussed is the interplay between, on 
the one hand, the intrinsic length and time scales of a soil and, on the 
other hand, the length and time scales associated with geometry and 
boundary conditions or emerging as consequences of a flow process. 
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